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Abstract 

We have combined Rossi X-ray Timing Explorer observations of X2127+119, the 
low-mass X-ray binary in the globular cluster M15, with archival X-ray lightcurves 
to study the stability of the 17.1 hr orbital period. We find that the data cannot 
be fit by the Ilovaisky et al. (1993) ephemeris, and requires either a 7a change 
to the period or a period derivative P/P ~ 9 x 10^'yr^ 1 . Given its remarkably 
low Lx / L opt such a P lends support to models that require super-Eddington mass 
transfer in a q ~ 1 binary. 
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1 Introduction 


The X-ray source X2127+119 is one of the ~10 bright (> 10 36 ergs _1 ) low 
mass X-ray binaries (LMXBs) located within globular clusters, being within 
2" of the core of M15. It was the first cluster X-ray source to have an optical 
counterpart proposed, the V~15 star AC211 (Auriere et al. 1984), on the ba¬ 
sis of blue colours, variability and location within the Einstein HRI error box 
(Geffert et al. 1989). The identification was spectroscopically conhrmed by 
Charles et al. (1986), who found He n A4686 and Ha emission. AC211 is opti¬ 
cally one of the brightest of all LMXBs, which is remarkable given its relatively 
low Lx- Consequently, it has been extensively studied (Hertz 1987, Callanan 

1988, Dotani et al. 1990, Ilovaisky et al. 1987, Naylor et al. 1988, Bailyn et al. 

1989, Ilovaisky et al. 1993) both in X-rays and at optical/UV wavelengths. 
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However, deriving the system parameters and understanding the geometry 
has proven difficult. Observations have yielded a wide variety of temporal be¬ 
haviour, obscuring the true binary period, and other unusual phenomena, as 
detailed below. 

AC211 exhibits the largest amplitude optical variations of any LMXB, and 
~5 times greater than in X-rays. Following their first photometric campaign 
Ilovaisky et al. (1987) detected a 8.54hr periodic modulation in their series 
of U-band CCD images, although there were puzzling unmodulated episodes. 
A reanalysis of HEAO-1 scanning data from November 1977 by Hertz (1987) 
revealed a consistent 8 . 66 hr X-ray modulation. However, further optical spec¬ 
troscopy by Naylor et al. (1988) gave a slightly different binary period of 
9.1 ± 0.5hr based on a radial velocity study of the He i A4471 absorption 
line. More remarkably, they found that the 7 velocity was blue-shifted by 
~150 ± lOkms -1 with respect to the cluster. Their favoured model for the 
system consisted of a high-inclination accretion disc corona source, which ac¬ 
counts for the unusually low Lx/L opt ~ 20 . A variable height disc rim (greatest 
when close to the accretion stream impact point) would then modulate both 
the X-rays and optical emission from the X-ray irradiated disc. They invoked 
ejection of the system from the cluster core (following a close encounter with 
another star) to explain the high 7 velocity. Alternative explanations were 
suggested; Fabian et al. (1987) proposed He 1 absorption in a supersonic wind 
from a massive corona, whilst Bailyn et al. (1989) preferred an outflow from 
the outer Lagrangian point as the site, requiring a common envelope system 
(the result of unstable mass transfer). The detection of a radius expansion 
X-ray burst with Ginga (Dotani et al. 1990) further complicated the picture, 
as clearly, in spite of the corona, the X-ray source must be directly visible 
for at least part of the time. Notably, this remains one of the most luminous 
X-ray bursts ever detected. 

The most recent piece in the puzzle arose from the continuation of Ilovaisky 
and co-workers’ photometric campaign. Ilovaisky et al. (1993) (hereafter 1193) 
analysed all their U-band data from 1984-89 as a single set, using a phase dis¬ 
persion minimisation (PDM) search method (Stcllingwerf 1978). This yielded 
a preferred period of 17.11 hr, double that previously accepted. This was able 
to explain the strange unmodulated observations and was found to be con¬ 
sistent with all the X-ray data up to that date. Using the minima in their 
lightcurve combined with that of the 1977 HEAO-1 data gave an ephemeris 
of <j) — 0.0 for JD 2447790.963 ± 0.018 + n x (0.713014d ± O.OOOOOld). 

Nevertheless, the AC 211 system remains an enigma. A high mass transfer 
rate, leading to significant mass loss from the system is a distinct possibility. 
Important constraints can therefore be placed by the measurement of any 
corresponding period change. Now that 9 years have elapsed since the last 
period and ephemeris determination we decided to make use of the opportunity 
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Table 1 

A Journal of the X-ray/optical observations of X2127+119/AC211. 


Observatory 

Detector 

Date (UT) 

Start End 

Time 

span 

Energy 

Range (keV) 

L X ’ ( 10kpc) 
(10 36 ergs -1 ) 

HEAO-1 

A-l 

19.11.77 

23.11.77 

96.4 lu¬ 

0.25-25 


~9.5 

EXOSAT 

ME 

30.06.84 

30.06.88 

ll.48 hr 

0.9-9 


~8 



22.10.84 

22.10.84 

15.27 hr 



~8 



20.10.85 

20.10.85 

7.17 hr 



~8 

Ginga 

LAC 

20.10.88 

24.10.88 

84.0 hr 

1-28 


6.4 

ASCA 

SIS 

16.05.95 

17.05.95 

26.3 hr 

0.4-12 


4.2 

ASCA 

GIS 

16.05.95 

17.05.95 

26.3 hr 

0.7-10 


4.2 

RXTE 

ASM 

23.02.96 

14.03.97 

384.6 d 

2-10 


4.2 



14.03.97 

02.04.98 

384.5 d 



3.4 

Optical 

CCD 

28.07.84 

22.09.89 

1881.4 d 

U-band 




presented by over 2 years of Rossi X-ray Timing Explorer All SKY Monitor 
(RXTE/ ASM) observations to revisit AC211. 


2 Determination of the revised ephemeris 


2.1 Data acquisition 


The HEASARC maintains comprehensive databases of past and present X- 
ray missions, from which we obtained the EXOSAT , AS'C'A/SIS and GIS 
lightcurves and the definitive dwell RXTE/ ASM data on X2127+119. Un¬ 
fortunately, the HEAO-1 and Ginga products are not readily usable, hence we 
simply digitised the data presented in the published lightcurve plots of Hertz 
(1987) and van Paradijs et al. (1990) (their Fig. 1 in both cases). Lastly, the 
CCD U-band data were only available in their published form, so we digitised 
the folded data as presented in Fig.8 of 1193. In spite of the relatively poor 
temporal sampling of the RXTE/ ASM dwell data, we found that the extensive 
coverage provided good phase sampling if taken over a sufficiently long time 
period. We chose to consider these data in two halves (as divided by time), 
hereafter referred to as ASM1 and ASM2. A log of all the data used is given 
in Table 1. 
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Fig. 1. Upper and lowest left panels: X-ray and optical lightcurves of 
X2127+119/AC211 respectively, folded (and binned) on the 1193 ephemeris and 
period of 0.713014 d. Lowest right panel: Sample of the X-ray lightcurves overlaid 
in chronological order (earliest lowest), with the flux units normalised and a vertical 
offset of 0.26 introduced between each dataset for clarity. Note the clear change of 
minimum phase over time. 

2.2 Measurement of the ephemerides 


All the data were folded on the ephemeris and period given in 1193, and binned, 
in order to increase the signal-to-noise. We used 40 phase bins, except for the 
HEAO-1 data (where we used 20), with each bin covering 0.05 (or 0.1) in 
phase and hence each overlapping the adjacent bins by ±0.025 (±0.05). This 
provided both good phase resolution, and some smoothing to reduce noise. 
These folded lightcurves are presented in Fig. 1. 


As expected, given its probably different physical origins, the optical lightcurve 
(labelled CCD) is different in shape to the X-ray, although both clearly demon- 
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Fig. 2. 0-C plot for the orbital ephemerides of X2127+119/AC211 for each X-ray 
observation and the one optical dataset; as determined by a Gaussian fit to each 
minimum (upper panel), and cross-correlations (lower panel). 


strate the eclipse by the secondary star. However, it is clear that the X-ray 
lightcurves also show considerable morphological evolution over time. The en¬ 
ergy ranges (see Table 1) are in fact similar for all the X-ray observations and 
so the changes must be intrinsic to the source. Whilst the EXOSAT and ASM 
data are similar, the intervening Ginga observation shows a totally different 
shape. 

In order to determine the relative orbital ephemerides for each observation, 
we adopted the same approach as 1193. They found the extended eclipse, or 
minimum, to be the only stable feature of the lightcurve, and so we used a 
simple Gaussian fit to determine the time of minimum in each case. Due to 
the poor phase coverage of the minimum in the EXOSAT data, we were forced 
to base our estimate on the sharp egress from minimum flux. The results are 
plotted in the form of an O-C plot (i.e. the phase offset between the expected 
minimum and that observed, versus cycle number) in the upper panel of Fig. 2. 
As a check, we also performed a cross correlation of the EXOSAT , ASM1 and 
ASM2 folds with the complete ASM lightcurve fold to determine their relative 
phase offsets (see lower panel of Fig.2; the values have been offset, such that 
the mean ASM results agree with those of the upper panel). 
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2.3 Evidence for a period derivative 


A variety of model fits have been made to the O-C data: (i) a constant O-C 
corresponding to a linear ephcmeris with the 1193 period of 0.713014 d is a 
very poor fit with xl — 27.1; (ii) a linearly increasing O-C corresponding to a 
linear ephemeris but with a different period of 0.713021 d (7 a larger then 1193) 
is statistically a good fit with xl — 1-12; (iii) a parabolic O-C corresponding 
to an ephemeris with the 1193 period and a period derivative of P = 1.75 ± 
0.90 x 10 -9 is also a good fit with xl — 0.59. However, one may apply a one¬ 
sided F-test (Bevington & Robinson 1992), to test the improvement of fit for 
an additional polynomial term between models (ii) and (iii). F x = Ax 2 /xl — 
(6.752 — 2.954)/0.59 = 6.43 corresponding to a 95% confidence level that the 
parabolic fit is better than the linear one. Moreover, the hypothesis of a linear 
ephemeris is undermined by the need for a 7a increase in the period, which 
is apparent in the discrepancy between the model and the HEAO-1 point. 
Comparison with the cross correlation results (lower panel) confirms that the 
data require a significant change to the 1193 ephemeris, whilst demonstrating 
that the ephemeris change derived from the minima times is if anything on 
the conservative side. We therefore consider that the evidence points towards 
a constant period derivative, with P/P ~ 9 x 10 _7 yr _1 . 


3 Discussion 


3.1 Physical origins of the X-ray and optical modulations 


It is profitable to compare the lightcurve of X2127+119/AC211 with those of 
the dipping sources X1822-371 and X0748-676 (Parmar et al. 1986), as shown 
in Fig.3. The X-ray lightcurve of X1822-371 clearly shows a broad minimum 
preceding, by 0.25 in phase, a deeper partial eclipse by the companion star, 
which is coincident with the optical minimum. The eclipse is only partial as 
the central source of this high inclination system is hidden by the disc and the 
X-ray scattering ADC is never fully occulted. To model the broad dip, White 
& Holt (1982) required a thickened region of the disc ~ 100° upstream of the 
line joining the companion and compact object. X0748-676 is essentially very 
similar, but has a lower inclination, thereby enabling the central source to 
be seen. Here, the almost total X-ray eclipse is always preceded by complex 
dip structure, which varies from cycle to cycle, although there appears to 
be a stable dip present at about </> = 0.65. This latter feature is once again 
probably due to a thickened region of the disc. Considering the lightcurves of 
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Fig. 3. The folded 1-10 keV EXOSAT /ME lightcurves of X1822-371 and X0748-676, 
using the period and ephemerides as given in White et al. (1981) and Parmar et al. 
(1986) respectively. The count rates have been normalized and the curves offset by 
1 for clarity. 


X2127+119/AC211 in Fig.l, the primary minima of every X-ray dataset are 
coincident with the optical minimum, but intervals of lower flux of variable 
duration are visible during phases </> = 0.4 — 1.0. This varying dip structure is 
certainly similar to X0748-676. The origin of the raised rim of the accretion 
disc is most probably the impact of the accretion stream. The understanding 
of this process remains incomplete, with a number of possible results from 
simulations, as outlined in Armitage & Livio (1998) and references therein. In 
addition to a bulge at (j) = 0.8, where the initial impact occurs, a reimpact 
following deflection and then reconvergence of the stream at (f) — 0.6 has been 
seen (Lubow 1989), but the azimuthal extent of these regions is unclear. 

If this interpretation is correct, then care clearly must be taken when deter¬ 
mining the ephemeris of the binary motion from the times of primary minima, 
as the dip structure will contribute to a varying degree to the flux reduction. 
This may well account for the changes seen in the asymmetry of the minima. 
Further observations, in order to both increase the base line of the ephemeris 
and to better model the details of the X-ray modulation are obviously needed 
to more accurately measure the period derivative. 
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3.2 Inferred mass transfer rate 


For this last section we will assume that P is within at least an order of 
magnitude of our estimate. The proposed evolutionary scenario for AC211 
involves a neutron star capturing a star close to the main sequence turn-off 
point in M15, which soon evolves towards the giant branch and turns on mass 
transfer by overflowing its Roche lobe (Bailyn & Grindlay 1987). Hence, one 
may apply the properties of a stripped giant to the secondary, which depend 
only on the mass of its helium core M c and not its total mass M 2 . 


The secondary radius R 2 (in R & ) and luminosity L 2 (in L & ) are given by the 
relations (King 1988): 


R'2 = 12.55 
L 2 = 33 ( 


(Me V 5 - 1 
VO.25/ ’ 

M c \ 8 ' u 
0257 


( 1 ) 

( 2 ) 


where M c is in M 0 and similarly for M 2 and Mx below. Using the Paczynski 
(1971) approximation for the Roche lobe radius: 


Rl 


0.462a 


M 2 


M 2 + Mx 


1/3 


(3) 


then eliminating the binary separation a by use of Kepler’s law and applying 
R 2 = Rl for the mass transfer condition yields: 


Mr 


7.65 / p s 

0.257 0 '°609 (d) 


(4) 


Taking the time derivative and dividing by (4): 


Mo 


.Mr 


P 


0-5—-r + 7.65—^ = - 


M, 


Mr 


P 


(5) 


The rate of change of M c with respect to M 2 is governed by nuclear evolution 
and for q = M 2 /M x < | (King 1988): 

-Mi = (6 ) 

M 2 | - 2 qM c y J 


Hence, we need to estimate the mass ratio q. There are two limiting cases 
for M 2 , when M c = M 2 and M c = 0.17 M 2 . The former is obvious whilst 



the latter follows from the need for M c /M 2 0.17, which is the Schonberg- 
Chandrasekhar limit, if the secondary is to have left the main sequence. Substi¬ 
tution into equation (4), with the 0.713014 d period gives M 2min = 0.146M 0 
and M'Z max = 0.92 M 0 . Thus the range of masses given in 1193 of 0.74—0.81 M 0 
(based upon the M15 isochrones and the appropriate age) fit well within these 
limits, towards the upper end. 


Finally, taking M 2 = 0.8 M 0 as a representative value and a canonical 1.4 M 0 
neutron star (q = 0.57) yields: 


M 2 _ _ 7 m c 
m 2 - 7 m c ’ 


(7) 


M 2 = -0.77 M 2 - 


( 8 ) 


and a mass transfer rate M 2 ~ — 6 x 10 7 Mq yr 1 . 

The corresponding accretion luminosity is given by: 

r (GM x M 2 \ Iif7 in46 ^ 

Lacc = -V -H- = — r 7 l.IT X 10 1 


Mo 


ergs 


-1 


v R x ) '' V M ©y r x j 

77 . 6 x 10 39 ergs _1 for M 2 ~ -6 x 10 -7 M 0 yr _1 


(9) 


The mass transfer rate is clearly super Eddington (L E ~ 77 .10 38 ergs _1 for a 
neutron star accreting hydrogen) and hence the accretion efficiency 77 -C 1. 
One would then expect significant mass loss from the system, after formation 
of a common envelope. This is just the scenario envisaged by Bailyn et al. 
(1989), and indeed our estimate for M 2 above is comparable to theirs, which 
was based on the He 1 absorption features. 


4 Conclusions 


We have obtained X-ray and optical lightcurves of X2127+119/AC211, from 
the HEASARC archives and published data, which span a total of 21 years, 
from 1977 HEAO-1 data up to the most recent 1998 RXTE/ASM data. Fold¬ 
ing the lightcurves on the 1193 ephcmeris and 17.1 hr orbital period, we have 
examined the stability of this period. The lightcurves in Fig. 1 clearly demon¬ 
strate the variety of X-ray morphologies, most probably arising from changes 
in the vertical structure of the outer edge of the accretion disc, which is re¬ 
sponsible for obscuring the X-ray emitting central region. As a result caution 
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is needed when using the folded lightcurves to check for any phase changes. 
However, the principal minimum (due to occupation by the secondary star) is 
obviously a persistent feature, and we have followed the procedure of 1193, to 
measure the relative ephemerides of each curve, based on its phasing. Only a 
non-zero period derivative applied to the original period, provides a good fit 
to these data, implying P/P ~ 9 x 10 -7 yr -1 . Modelling the secondary as an 
evolved sub-giant, it is possible to estimate the mass transfer rate necessary 
to drive such a period change. This yields M 2 ~ — 6 x 10~ 7 M 0 yr^ 1 , and hence 
an accretion luminosity L acc ~ rj. 6 x 10 39 ergs -1 . Lastly, as this far exceeds 
the Eddington luminosity, the accretion efficiency must be very low , requiring 
there to be significant mass loss from the system. This is consistent with the 
earlier spectroscopic data of Bailyn et al. (1989), and a q ~ 1 binary. 


Acknowledgement 


This research has made use of data obtained through the High Energy As¬ 
trophysics Science Archive Research Center Online Service, provided by the 
NASA/Goddard Space Flight Center. In particular, we thank the ASM/RXTE 
teams at MIT and at the RXTE SOF and GOF at NASA’s GSFC for provision 
of the ASM data. LH acknowledges support of a PPARC studentship. 


References 

Armitage P J , Livio M 1998 , ApJ. , 493, 898. 

Auriere M, le Fevre O , Terzan A 1984 , A&A. , 138, 415. 

Bailyn C D, Garcia M R , Grindlay J E 1989 , ApJ. , 344, 786. 

Bailyn C D , Grindlay J E 1987 , ApJL. , 316, 25. 

Bevington P R , Robinson D B 1992 , Data Reduction and Error Analysis for the 
Physical Sciences McGraw-Hill. 

Callanan P J 1988 , Ir. Astr. J. , 18, 152. 

Charles P A, Jones D C , Naylor T 1986 , Nature , 323, 417. 

Dotani T et al. 1990 , Nature , 347, 534. 

Fabian A C, Guilbert P W , Callanan P J 1987 , MNRAS. , 225, 29P. 

Geffert M, Auriere M, Ilovaisky S A , Terzan A 1989 , A&A. , 209, 423. 

Hertz P 1987 , ApJL. , 315, 119. 


10 



Ilovaisky S A, Auriere M, Koch-Miramond L, Chevalier C , Cordoni J P 1987 , 
A&A. , 179, LI. 

Ilovaisky S A, Auriere M, Koch-Miramond L, Chevalier C, Cordoni J P , Crowe 
R A 1993 , A&A. , 270, 139. 

King A R 1988 , QJRAS. , 29, 1. 

Lubow S H 1989 , ApJ. , 340, 1064. 

Naylor T, Charles P A, Drew J E , Hassall B J M 1988 , MNRAS. , 233, 285. 
Paczynski B E 1971 , ARA&A. , 9, 183. 

Parrnar A, White N E, Giommi P , Gottwald M 1986 , ApJ. , 308, 199. 
Stellingwerf R F 1978 , ApJ. , 224, 953. 

van Paradijs J, Dotani T, Tanaka Y , Tsuru T 1990 , PASJ. , 42, 633. 

White N E , Holt S S 1982 , ApJ. , 257, 318. 


11 



